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Abstract 

We identify a gravitational-dynamical process in near-Keplerian potentials of galactic nuclei that occurs 
when an intermediate-mass black hole (IMBH) is migrating on an eccentric orbit through the stellar cluster 
towards the central supermassive black hole. We find that, apart from conventional dynamical friction, 
the IMBH experiences an often much stronger systematic torque due to the secular (i.e., orbit-averaged) 
interactions with the cluster's stars. The force which results in this torque is applied, counterintuitively, in 
the same direction as the IMBH's precession and we refer to its action as "secular dynamical anti-friction" 
(SDAF). We argue that SDAF, and not the gravitational ejection of stars, is responsible for the IMBH's 
eccentricity increase seen in the initial stages of previous A-body simulations. Our numerical experiments, 
supported by qualitative arguments, demonstrate that (1) when the IMBH's precession direction is artificially 
reversed, the torque changes sign as well, which decreases the orbital eccentricity; (2) the rate of eccentricity 
growth is sensitive to the IMBH migration rate, with zero systematic eccentricity growth for an IMBH whose 
orbit is artificially prevented from inward migration; and (3) SDAF is the strongest when the central star cluster 
is rapidly rotating. This leads to eccentricity growth/decrease for the clusters rotating in the opposite/same 
direction relative to the IMBH's orbital motion. 

Subject headings: black hole physics — Galaxy: center — stars: kinematics and dynamics 



1. INTRODUCTION 

Recent studies on secular dynami cs in near-Keplerian po- 
tentials dRauch & Tremairiel 119961: iHopman & Alexander! 

it 



20061: lEilon etalj l2009t iKocsis & Tremaind 1201 L 
Madig an et alJl20l! ) have focused attention on the long-term 



dynamical evolution of stars and compact objects in galactic 
nuclei. Due to persistent gravitational torques exerted by 
stellar orbits on each other as they precess slowly around the 
supermassive black hole (SMBH), the angular momentum 
evolution of stellar orbits can proceed at a much higher 
rate than that of energ y evolution (resonant relaxation; 
iRauch & Tremainelll996l) . In this paper we extend this re- 
search to study secular dynamical effects on the inspiral of a 
massive body, e.g., an intermediate-mass black hole (IMBH), 
into the combined potential of an SMBH and nuclear stellar 
cluster. 

As a body of mass M moves through a distribution of field 
stars of individual mass m a , it experiences a frictional force 
anti-parallel to its velocity v, 



M— 

dt 



-16ir 2 G 2 M 2 m a In A 



dv a vlf{v a ) 



v 

(1) 



This is Chandrasekh a r's dynamical frictio n f ormul a 

dChandrasekhari 119431: iTremaine & Weinbergl 11984b for 
an isotropic field star distribution normalized such that the 
density n(r) = J d 3 v a f(r,v a ), and In A is the Coulomb 
logarithm. The formula accounts only for stars moving 
slower than the massive body, and it neglect s self-g ravity 
of the stars themselves; lAntonini & M erritt (2012) have 
recently updated the formula to include changes induced 
by the massive body on the stellar velocity distribution and 
the contribution from stars moving faster. The stars are 
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deflected by the massive body into a gravitational "wake" 
behind it which results in an increase in st ellar density, the 
amplitude of which is proportiona l to M (iDanbv & Cam ml 
119571; iKalnai si [1971 lMulderJ[T983l) : hence, the gravitational 
force experienced by the massive body is proportional to M 2 . 
This frictional force acts to decrease the kinetic energy and 
angular momentum of the massive body and it sinks towards 
the center of the potential. 

The picture described above does not take into account the 
orbit-averaged dynamics that is particular for near-Keplerian 
potentials and that has been shown to play a central role in 
the angular momentum evolution of stars near the SMBH. In- 
deed, it is natural to assume that some other form of dynami- 
cal friction could be associated with the orbit-averaged poten- 
tial created by the IMBH, as this potential is rotating around 
the SMBH due to precession of the IMBH's eccentric orbit. 
The possibilit y of this secular dyn amical friction was already 
suggested by IRauch & Tremaind in 1996. In this paper we 
explore this effect and find that, contrary to our experience 
with ordinary dynamical friction, the resulting torque comes 
from a gravitational force acting in the same direction as that 
of the IMBH precessionQ We therefore refer to this effect as 
"secular dynamical anti-friction" (SDAF). 

For the remainder of the paper we will keep referring to 
the massive in-spiraling body as an IMBH, though the reader 
should keep in mind that the essential dynamics will hold for 
any massive body moving in a near-Keplerian potential. As 
large mass ratio massive black hole binaries are expected to 
coalesce in at least some merging and non-merg ing galaxies 
jMilosavlievic & Merrittir2 003; Pre toet al.lfeoill) . the dynam- 
ics of such systems are important to understand. 

1 This is true if the stellar cluster is non-rotating. As we will see below, 
if the cluster is rapidly rotating, then the torquing gravitational force can be 
applied in the same direction as the cluster's rotation. 
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In particular, IMBH inspirals into SMBHs will be a ma- 
jor source of gravitational waves for space-based laser in- 
terferometers such as the proposed European New G ravita- 
tional Wave Observatory (lAmaro-Seoane etaD 120121) . Sim- 
ulations of the inspiral of IMBHs (10 2 - 1O 4 M ) through 
nuclear stellar clusters embedding an SM BH have shown 
an in cr ease in eccentricity of the IMB H (Baumgar dt et al.l 
20061: Pat subavash ietaLl 120071; lLockmann et al.1 120081: 
Iwasawa et alj 1201 U [Sesana et al.l 1201 U lAntonini & M erritt 
20121: iMeiron & Laorll201 ll) . Here, we show that the theory 



of SDAF can explain the initial stages of this phenomenon 
and present results of A-body simulations set up to test this 
claim. 

2. SECULAR DYNAMICAL ANTI-FRICTION 

Let us consider the combined gravitational potential of an 
SMBH, mass M., embedded in a nuclear star cluster with 
a power-law number density distribution, n(r) oc r~ a , and 
individual stellar masses m. The specific energy and angular 
momentum of a stellar orbit can be written in terms of Kepler 
elements as 

B-™ (2, 



2a 



| J | = |r x v| = [GM.a(l - e 2 )] 



2N11/2 



(3) 



where a is the semi-major axis of the elliptical orbit, and e 
is the eccentricity. The eccentricity vector, e, of the star is 
that which points from the occupied focus of the orbit to the 
periapsis of the orbit and has a magnitude equal to the scalar 
eccentricity of the orbit, 



1 



GM. 



(v x J) 



(4) 



The orbit of a star in this potential is nearly Keplerian (i.e., 
closed) but e will precess due to the additional (Newtonian) 
potential and general relativity. We define the precession time, 
iprcc, as a timescale over which an orbit precesses by 2tt rad 
in its plane. As Newtonian precession acts with retrograde 
motion (i.e., in the direction opposite of a star's motion) and 
general relativistic precession with prograde motion, 



1 



1 



fd +GR 
prcc prcc 



The general relativistic precession time is given by 
o(l 



t 



GR 
prec 



e 2 )c 2 



3GM. 



P(a), 



(5) 



(6) 



with the orbital period P(a) = 2ny / a 3 /GM m . The New- 
tonian precession time due to the addition of the star cluster 
mass in an otherwise-Keplerian potential, tp rec , can be written 

for a 7^ 2 as 



=7r(2-a)/(e,a)- 1 



M. 

N7m 



P(a) 



(7) 



where A< is the number of stars within semi-major axis a, 
and 



/(6 ' a) = ^ [jj tl 



J 2 



l + \fT~J 2 
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FIG. 1. — Dependence of the precession time (M. / (N < m)P(a) = 1) 
of a stellar orbit on orbital eccentricity, e. Several values are plotted for 
different power-law densi ty profiles, n(r) oc r~ a . The dependence given in 
IMerritt & Vasiliev (2011, Mil) and llvanov etal] j200l 105) are plotted for 
comparison. 



(lLandau & Lifshitzl fl969h . Here J = J/fGM.a) 1 / 2 = 
(1 — e 2 ) 1 / 2 . For all values of a, the Newtonian precession 
time of a stellar orbit increases with eccentricity. The dif- 
ference in precession times for circular and highly eccentric 
stellar orbits varies with a. We plot t prcc in Figure Q] set- 
ting the quantity in brackets [M./(A <-m)P(a)1 in Equatio n 
(0 equal to 1, with /(e, a) as given in|Madigan et al.l d201 ll) . 
For a steep cusp, a = 1.75 dBahcall & Wold 119761) . there 
is a factor of ~ 5 difference in precession time for a near- 
circular (e = 0.01) and near-radial (e = 0.99) orbit, which 
increases to ~ 7 for a shallow profi le a = 0.5. We plot Equa- 
tion (13) from IMerritt et al.l (120111) and Equation (Al 1) from 
llvanov et all (120051) (a = 1.5) as a comparison. Differences 
arise in the normalization of the functions due to approxima- 
tions made in deriving the analytical formulae in these papers; 
the monotonic dependence of the precession time on orbital 
eccentricity persists. We note that the e xact expressions for 
the precession time in llvanov et al.l (120051) agrees with numer- 
ical evaluation of our Equation dT). 

We now introduce an IMBH of mass Mj m bh on an orbit 
with non-zero eccentricity, ei m bh, and semi-major axis, ai m t>h, 
within the dynamical radius, r^, of the SMBH, defined such 
that the mass in stars at equals that of the SMBH, iV"(< 
r/j)m = Al,. As we are interested in secular dynamics^ it 
is useful to envisage the mass of the IMBH spread smoothly 
out along its orbit, such that the local density at any segment 
is inversely proportional to its local velocity. The position of 
furthest distance from the SMBH, the apoapsis of the orbit, 
will therefore contain the most mass. Over a timescale longer 
than its orbital period but much less than the precession time, 
it will exert a strong (specific) gravitation al torque on the orbit 
of a s tar with a similar semi-major axis (iGiirkan & Hopmanl 
2007) on the order of 



M = |r xFl 



GMimbheimbhe , 

Oqmbh 



(9) 



(8) 



where dip is the angle between two orbits. 

2 The term "secular" refers to long-period dynamics, in which the depen- 
dence on the mean longitude of an orbi t is dropped from the disturbing func- 
tion (see, e.g., Murray & Dermott 1999). 
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FIG. 2. — Left: schematic of orbits of IMBH (black, filled line) and prograde leading (orange, dashed line) and prograde trailing (blue, dotted line) stellar orbits 
in the frame of reference rotating with the IMBH orbit. The velocity vectors of each are indicated at apoapsis and the direction of precession of the IMBH is 
shown with a curved arrow at top left. Right: schematic of orbits of IMBH (black, filled line) and retrograde leading (orange, dashed line) and retrograde trailing 
stellar orbits (blue, dotted line) in the frame of reference rotating with the IMBH orbit. 



The orbit of the IMBH will not affect the energy (or equiv- 
alently, the semi-major axis) of the stellar orbit as its gravita- 
tional potential is stationary over this timescale. A star on a 
coplanar orbit with the IMBH will experience a torque that is 
parallel (or anti-parallel) to its angular momentum vector J. 

Hence, the torque will act to change the magnitude of J; 
equivalently, the eccentricity of the orbit. Let us consider the 
stellar orbit with an eccentricity e ~ ej m bh and semi-major 
axis a ~ aj m bh- We place the orbit such that its eccentricity 
vector forms an angle 5(f> with that of the IMBH. 

If there existed no gravitational interaction (secular or oth- 
erwise) between the IMBH and the star, they would precess 
in the background stellar potential at the same rate. Their ec- 
centricity vectors would rotate by 2ir rad in one i prcc and 5(f) 
would remain constant. Switching on secular gravitational in- 
teractions results in a strong secular encounter between the 
star and the IMBH. The star's orbit is pushed away from the 
IMBH, transferring negative or positive angular momentum 
to the orbit of the IMBH, depending on its orientation. We 
elaborate on the dynamics in the following. 

If the star is moving on a prograde orbit (in the direction of 
motion of the IMBH) that leads in front of the IMBH (in terms 
of precession direction; see Figure [2] left), it will experience 
a torque in the direction of its angular momentum vector J 
such that dJ/dt > 0. Due to the resulting decrease in eccen- 
tricity, the stellar orbit will precess faster (see Figure [1} and 
move away from the IMBHQ Conversely if the stellar orbit is 
trailing, the torque exerted on it will be anti -parallel to J. The 
stellar orbit will increase in eccentricity, precess more slowly 
and in effect be pushed away from the orbit of the IMBH. In 
both cases, d\S(f>\/dt > 0. 

A star moving along a retrograde orbit (see Figure |2] right) 
precesses in the opposite direction to that of the IMBH. If 
the orbit is leading in front of the IMBH, the torque will act 
to increase its eccentricity and hence decrease the differen- 
tial precession rate between the two orbits. This increases the 
timescale over which the stellar orbit experiences a coherent 
torque and its eccentricity can increase to one. In this scenario 
the orbit can easily "flip" through e = 1 to a prograde orbit, 

3 This mechanism involves the same inteiplay between stellar torques and 
orbital precession which r esults in an instability in an eccentric disk of stars 
near an SMBH I Madigan et aU200fJD , 



and be pushed away from the IMBH as described above. If if 
does not flip, the retrograde orbit will precess past the IMBH 
orbit, experience a torque which will decrease its eccentric- 
ity and hence will increase the differential precession rate be- 
tween the two orbits. Once again, in both cases the stellar 
orbit is pushed away from that of the IMBH, d\6<f>\/dt > 0. 

In contrast to Chandrasekhar's dynamical friction, which 
acts to increase the density of stars behind a massive body, 
the effect of SDAF is to push stellar orbits away from that of 
the massive body. The varying strength of Chandrasekhar's 
dynamical friction on different segments of a massive body's 
orbit causes evolution of orbital eccentricity, the sign of which 
is dependent on the background stellar d ensity profile (e.g., 
iGould & Ouilienll200l iLevin et alj|2005l) . With SDAF, it is 
the relative number of prograde/retrograde stellar orbits which 
determines the sign of the orbital eccentricity evolution of the 
IMBH. 



2.1. Increase in orbital eccentricity 

During a strong secular encounter of the IMBH with a star, 
the angular momentum of the star is randomized. The rate of 
encounters of the IMBH with retrograde stars is higher than 
that of prograde stars, and the net result is to decrease the 
angular momentum of the IMBH. If the precession direction 
of the IMBH is artificially reversed, the net torque changes 
sign and the angular momentum of the IMBH increases. Mi- 
gration is important in this picture. If the IMBH does not 
migrate, then the relative number of prograde and retrograde 
orbits quickly adjusts so there is no net secular torque in the 
steady state. However, if the IMBH migrates due to dynam- 
ical friction, the adjustment does not have time to occur and 
there is an anti-friction torque. In the next section, we perform 
numerical simulations to test this picture. 

2.2. Comparison with theories in the literature 

The iV-body simulations of Iwas awa et al.l (1201 lh show an 
increase in the eccentricity of an IMBH as it spirals into 
an SMBH. With careful analysis of the orbital properties of 
the field stars, the authors attribute the increase to secular 
chaos (and the Kozai mechanism), brought about by the non- 
axisymmetric potential induced by the IMBH. 
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This acts to bring the relative number of orbits back into 
balance after prograde-moving stars are preferentially ej ected 
by the IMBH (retrograde stars have larger relative velocities 
with respect to the IMBH than prograde, and hence are ejected 
less frequently). The retrograde orbits moving to prograde 
orbits extract angular momentum from the orbit of the IMBH 
a nd its eccentricity in creases. 

ISesana et al.l d201 lb present both hybricQ and direct A -body 
simulations with stellar cusps of varying degrees of rotation, 
to examine the importance of the relative number of pro- 
grade/retrograde orbits on the orbital evolution of the IMBH. 
The authors propose a hypothesis in which the ejections of 
stars, via the slingshot mechanism, are the cause of the evolu- 
tion of the IMBH eccentricity. Most retrograde orbits which 
end up being ejected from the potential, donate their negative 
angular momentum to the IMBH and move to prograde orbits 
before they escape. Again it is the preferential ejection of pro- 
grade orbits that drives the increase of the orbital eccentricity 
of the IMBH. 

With SDAF, it is the cumulative secular torques from retro- 
grade stellar orbits that lead to the increase in orbital eccen- 
tricity of an IMBH, irrespective of ejections of stars from the 
potential. SDAF makes testable predictions which can distin- 
guish itself from the above theories. 

1 . If the orbit of the IMBH is made to artificially precess 
in the opposite direction to its true motion, its rate of 
change of angular momentum should reverse sign as 
well. 

2. If escaping stars are injected back into the cluster with 
the same parameters as before their interaction with the 
IMBH, the eccentricity evolution of the IMBH should 
not qualitatively be affected. 

3. Due to the "pushing away" of coplanar stellar orbits, 
there will be a density anisotropy in the stellar distribu- 
tion as a result of SDAF. 

4. If dynamical friction is artificially slowed, there will be 
a lower rate of axisymmetric flow of retrograde orbits to 
one side of the IMBH's orbit and hence a lower negative 
rate of change of angular momentum. 

In the following section we present the results of numeri- 
cal trials to test the theory of SDAF and decrease of orbital 
angular momentum of an IMBH in a galactic nucleus. 

3. NUMERICAL METHOD AND SIMULATIONS 

The simulations presented in this paper are performed using 
a special-purpose, mixed-variable symple ctic A-body inte- 
grator as described in Madigan et al.ld2Qllb . We use a galactic 
nucleus model with parameters chosen to represent the Galac- 
tic center. At the center of the coo rdinate system, th ere is an 
SMBH with M. = 4 x 1O 6 M0 dGhez et all 120031) . and an 
IMBH with a mass M imb h within the range [10 3 - 1O 4 ]M 
with 0.04 pc < a < 0.08 pc. We simulate the nuclear star 
cluster as a smooth potential with a power-law number den- 
sity profile n(r) oc r~ a ; this decreases the computation time 
and lends a clear, direct interpretation to our results. Typically 
a E [0.5,1.75], and is normalized such that at one parsec 

4 The hybrid model couples numerical three-body scatterings to an analyt- 
ical formalism for interactions between the cusp and the IMBH I Ses ana et alj 
2008). 
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FIG. 3. — Example of a secular horseshoe orbit of a star of mass m = 
1M Q , bound to an SMBH of mass M. = 4 X 1O 6 M . The stellar orbit 
feels a gravitational torque from an IMBH (mass M; m t,h = 10 Mq) and 
precesses due to general relativity and a background stellar potential. Top: 
the eccentricity vectors of the IMBH and the prograde-directed single star 
in the frame co-precessing with the eccentricity vector of the IMBH. Bottom: 
eccentricity evolution of the stellar orbit over 8 Myr during which it performs 
several secular horseshoe orbits. 



the stellar mass is M(< lpc) = 1 x 10 6 M ^ dSchodeletalJ 
120071: iTrippe et al.l2008tlSchodel et aI1l2009l) . To test our the- 
ory of SDAF, we include an isotropically distributed popula- 
tion of test stars with individual masses m within the range 
[10— 100] Mq which gravitationally interact only with the 
IMBH, not with each other. In doing so we isolate the effects 
of the IMBH on the stellar orbits. The orbits of the IMBH 
and the test stars precess according to Equation (|5), such that 
for each time step, dt, they rotate in their orbital plane by an 
angle \5<f>\ = 2iz(dt/t wcc ). 

As the IMBH moves in a smooth stellar potential, it expe- 
riences no dynamical friction from it. To ac count for this we 
use the dynamical friction formulae given by lJust et al.l (1201 ll) 
to a rtificially decrease the semi-major axis of the IMBH's or- 
bit. lJust et all d201 ll) use a self-consistent velocity distribu- 
tion function in place of the standard Maxwellian assumption 
and an improved formula for the Coulomb logarithm to arrive 
at dynamical friction timescale for a near-Keplerian potential, 
which can differ by a factor of three from the standard formul a 
in Equation (Q]). 

The radial evolution of the IMBH (point-like object) be- 
haves, for a ^ 3/2, a^| 



r 1 



t \ 2Q " 3 

tj3F J 

5 For a = 3/2 see similar equations in Just et al. < 20 11). 



(10) 
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where £df is the dynamical friction timescale calculated from 
the starting position of the massive body ro, 



P 



Mi 



2 7 r(3-a)xlnAM imbh Af< 



(11) 



X is the fraction of stars with velocity smaller than that of 
the IMBH (circular case), 



X(a) 



1 



V2B{Z/2,a-l/2)j 
and A is given by 

(3 + a) M. 



a-3/2 

I - y 1 du 



A = 



a(l + a) M imbh ' 



(12) 



(13) 



We implement Equations ( TTOb - ( TT3T > by adjusting the semi- 
major axis of the IMBH orbit at each time step with the fol- 
lowing transformation 



j 9= 1- 



iDF(n-i) 



(14) 



(15) 



Essentially we are re-scaling the stellar ellipse, without 
changing its eccentricity. In the next subsections we present 
three examples of dynamics resulting from SDAF. 

3.1. Secular horseshoe orbits 

Horseshoe orbits describe the motion of satellite bodies in 
the circular planar restric ted three-body problem (f or a com- 
prehensive overview see Murray & Dermottlll999l) which li- 
brate on paths encompassing the L3, LA, and Lb Lagrangian 
points. Such orbits map out horseshoe shapes in the frame co- 
rotating with two massive bod ies and have been observ ed in 
the Saturnian satellite system (Dermott & Murray 198 Tb and 
in the co-orbital Earth Asteroid 2002 A A29 ((Connors et al. 
12001 . 

Secular dynamical anti-friction results in secular horse- 
shoe orbitsQ Secular horseshoe orbits differ from those 
in our above description in that it is the eccentricity 
vector, as opposed to the guiding center of the orbit 
which maps out a horseshoe in the frame co-rotating 
with the IMB H (for close analog ies in planetary dynam- 
ics, see iMichtchenko & Malhotrai 120041 iPan & Saril l2004t 
Lith wick& Wul 120111) . We illustrate this behavior in Fig- 
ure |3 where we show the position of the eccentricity vec- 
tors of a IMBH and a star in the frame of reference co- 
precessing with the eccentricity vector of the IMBH. In this 
example, M. = 4x 10 6 M Q , M imb h = 10 3 A1©, m = 1M Q , 
a = 1, M(< lpc) = 1 x 1O 6 M , a imbh = a = 0.0703 pc, 
eimbh = 0.7, and e = 0.69. Both the IMBH and stellar orbits 
lie on the xy-plane with an initial angle of 5<fr = 7r/4 rad be- 
tween their eccentricity vectors. This behaviour is transitory 
due to non-secular gravitational interactions. The timescale 
over which the star in this simulation performs a horseshoe in 
eccentricity space is ~ 1.5 Myr or ~ 1700 orbital periods. 

6 Referred to as "windshield-wiper orbits" in Merritt & Vasiliev 1 201 ID in 
the context of non-axisymmetric star clusters. 
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FIG. 4. — Schematic of octants in eccentricity vector space (top). The unit 
eccentricity vector of the IMBH lies parallel with (1,1,1) in octant 1 (black 
dot). Fractional change in number of prograde (middle) and retrograde (bot- 
tom) stars in each octant in the frame rotating with the IMBH. Time is given 
in units of initial orbital period of the IMBH, P(a = 0.06 pc). The increase 
in the number of prograde stellar orbits in octant 8, and retrograde stellar 
orbits in octant 1 , is as predicted. 

3.2. Separation of potential into pro- and retro- grade orbits 

As illustrated in Figure [3] the eccentricity vectors of pro- 
grade orbits which have a small inclination angle with respect 
to an IMBH will tend to librate about a point it rad opposite 
to that of the IMBH - the secular analog of the Lagrangian 
L3 point. In contrast, the eccentricity vectors of retrograde 
orbits (those which do not flip to prograde orbits) will, in a 
time-averaged sense, cluster around that of the IMBH before 
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they are pushed away, as they precess more slowly due to the 
increase in their orbital eccentricity. We find that the presence 
of a massive body, such as an IMBH, surrounded by less mas- 
sive objects in a near-Keplerian potential tends to segregate 
the lighter objects into a region dominated by prograde orbits 
and one dominated by retrograde. 

This mechanism has the strongest effect for, but is not 
restricted to, co-planar orbits, as scalar angular momentum 
changes result in eccentricity changes whereas vector angular 
momentum relaxation re-orients the stellar orbital plane. 

In Figure H] we plot an example of this segregation. The 
top figure shows a schematic of eight octants in eccentricity 
vector space. We rotate the cluster such that the eccentric- 
ity vector of the IMBH lies parallel with (1,1,1) and follow 
the change in the number of massless stellar orbits with ec- 
centricity vectors in each octant. The middle plot shows the 
results for prograde orbits and the bottom plot for retrograde. 
There is a clear increase in the number of prograde stellar 
orbits in octant 8, opposite that in which the IMBH lies, as 
expected. Furthermore, there is an enhancement of retrograde 
stellar orbits in octant 1 as the theory predicts. 

3.3. Increase of orbital eccentricity of IMBH 

Here we examine two of the main testable predictions of 
SDAF: the eccentricity evolution of an IMBH with (1) arti- 
ficially reversed apsidal precession and (2) with a decreased 
migration rate from Chandrasekhar's dynamical friction. We 
present the results of short-term simulations performed to test 
these predictions, using ten realizations for each simulation. 
In the following, M. = 4 x 1O 6 M , M imbh = 10 4 M Q , 
aimbh = 0.06 pc, e im bh = 0.4, a = 1.5, and M(< lpc) 
I x 1O 6 M . We integrate the orbits of 1000 test stars in 
each simulation, distributed from 0.006 pc < a < 0. 12 pc, with 
m = 20M Q . 

In Figure [5] we show the evolution of orbital eccentricity 
of the IMBH as a function of time in units of the initial or- 
bital period, P(a = 0.06 pc) ~ 692 yr. Each line shows the 
mean eccentricity of ten simulations. C proc is a pre-factor 
placed in front of the apsidal precession angle, 5<fi, at each 
time step. The reference simulation with the correct preces- 
sion rate C plC c = 1-00 shows an increase in eccentricity as 
expected. Increasing the precession rate by a factor of 1.25 
does not however increase the rate of eccentricity growth, 
as the IMBH orbit has less time to exert torques on the sur- 
rounding stellar orbits. In the next two simulations we reverse 
the direction of apsidal precession of the IMBH, keeping the 
precession period roughly constant in relation to the first two 
simulations. The overall eccentricity growth of the IMBH is 
negative, as predicted. 

Next, we present the results of simulations wherein we 
change the migration rate of the IMBH (Figure [6j. SDAF 
predicts that lower migration rates lead to slower eccentric- 
ity growth as the IMBH interacts with retrograde orbits at a 
lower rate. In this figure, Ctdf is a multiplicative pre-factor 
in front of the dynamical friction timescale given in Equa- 
tion ( fTTT i. The semi-major axis of the IMBH orbit shrinks 
from 0.06 pc to ~ 0.024 pc over the simulation for Ctdf = 1. 
For Ctdf = 2, 5, 10 the eccentricity growth of the IMBH 
is suppressed (though still positive), exactly as SDAF pre- 
dicts. In the simulation in which we half the friction timescale 
(Ctdf = 0.5), and hence increase the migration rate by a fac- 
tor of two, the eccentricity growth rate does not increase. We 
attribute this to the shortening over the timescale over which 
the IMBH exerts persistent torques on the surrounding stellar 
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FIG. 5. — Mean eccentricity evolution (each curve is the average of ten 
simulation runs) of an IMBH experiencing four different rates of apsidal pre- 
cession. The first is the reference simulation with the correct precession rate 
corresponding to the smooth background potential (red solid line). In the 
second simulation, the IMBH has an increased precession rate by a factor of 
1.25 (orange dotted line). In the third simulation, the IMBH has an increased 
precession rate by the same factor but the IMBH is made to precess in a re- 
versed direction (green dotted line). In the fourth simulation, the IMBH has 
a precession period corresponding to the reference simulation but again in 
a reversed direction (blue dotted line). Time is given in units of the initial 
orbital period P(a = 0.06 pc). 
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FIG. 6. — Mean eccentricity evolution (each curve is the average of ten sim- 
ulation runs) of an IMBH with five different migration rates due to dynamical 
friction. SDAF predicts lower rates of eccentricity growth with increasing dy- 
namical friction timescales. The parameter C t df is a multiplic ative pre-factor 
in front of the dynamical friction timescale given in Equation (IT) . 



orbits. 

Finally, we run simulations wherein we soften the IMBH's 
gravity to such an extent that it switches off entirely for a 
star coming within 0.01 pc of the IMBH itself. We do so to 
confirm that secular dynamics are dominating the eccentric- 
ity evolution of the IMBH orbit; these conditions inhibit the 
three-body scattering s and ejections which are the basis for 
the the ories given in llwasawa et al.l ( 1201 11) and iSesana et al.l 
d20llb . We find it makes no significant difference for the an- 
gular momentum evolution of the IMBH even though it sig- 
nificantly reduces the number of stellar ejections (by a factor 
of- 5). 
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FIG. 7. — Mean eccentricity evolution (each curve is the average of ten sim- 
ulation runs) of an IMBH orbiting in counter- and co-rotating stellar cusps. 
The evolution is faster than in non-rotating stellar cusps. 



The problem of dynamical friction on a bar rotating 
through a spherical s tellar system has been studied by 
Trema ine & Weinberg! ([1984). The authors find that the 
torque is directed opposite to the direction of the bar's rota- 
tion; see, e.g., their Equation (67). This seems in contradic- 
tion with our finding that the precessing IMBH orbit experi- 
ences a torque directed in the same direc tion as its precession. 
We no te however, that Equation (67) of Trem aine & Weinberg! 
( 1984) is derived under the assumption that the bar's mass is 
small enough such that its effect on the orbits of the stars is 
weak and perturbation theory can be used. We, on the other 
hand, explore the regime where the relevant secular interac- 
tions between the IMBH and the stars are strong. A secular 
horse-shoe orbit is an example of such an interaction, an an- 
gular momentum flip is another one. 

3.4. Eccentricity evolution in a rotating stellar cluster 

For a non-rotating stellar cusp, the eccentricity increase in 
the IMBH orbit due to SDAF is relatively small (as seen in 
Figures [5] and [6}. However, net ro tation is ob served in most 
well-r esolved galactic nu clei (e.g., Seth 2010), including our 
own dGenzel et al.lll996h . In this case, there is a higher rate 
of strong secular encounters between either retrograde or pro- 
grade stellar orbits and that of the IMBH and, as a result, the 
rate of eccentricity evolution of the IMBH orbit due to SDAF 
is greatly enhanced. Specifically, SDAF predicts a strong in- 
crease in orbital eccentricity of the IMBH in counter-rotating 
clusters due to the increased number of retrograde stellar or- 
bits, while it predicts a decrease for co-rotating clusters due 
to the increase in prograde orbits. 

In Figure 7, we show the results of simulations with co- 
and counter-rotating stellar cusps to demonstrate the signifi- 
cant change in eccentricity of the IMBH above that achieved 



in non-rotating cusps. Using the same parameters as in our 
previous simulations, we see a mean increase (decrease) in 
eccentricity of e = 0.4 to e w 0.66 (0.32) in the counter- 
(co-) rotating stellar cusp scenario. Demonstrating a more 
rapid growth in eccentricity requires additional "grainy" mass 
in surrounding cusp stars (i.e., A-body/field particles as op- 
pose d to a smooth poten tial). We refer the reader to Figure 
2 of ISesana et al.l d201 ll) who simulate a similar potential to 
that in our model but use significantly more A-body parti- 
cles. They find a growth (decrease) in eccentricity of e = 0.5 
to e sa 1 (0.06) ove r a similar timescale. However, in con- 
trast to ISesana et all s interpretation, we attribute this growth 
(decrease) to SDAF. 

4. DISCUSSION 

We have presented a new secular gravitational-dynamical 
process in galactic nuclei which acts to decrease the orbital 
angular momentum of a massive body. There are two essen- 
tial components to the decrease. The first is the net flow of 
retrograde orbits past that of the massive body which extracts 
orbital angular momentum. The second is Chandrasekhar's 
dynamical friction force acting on the massive body which 
decreases its semi-major axis and provides a mechanism for 
net negative angular momentum change. Although we have 
told this story in the context of an IMBH inspiral into the cen- 
ter of a galaxy, the interactions described here will be valid 
for any massive body in a near-Keplerian potential. 

SDAF can account for the rise in eccentricity only in the 
initial stages of the inspiral of the IMBH, not during the 
stalling phase. The eccentricity d oes rise strongly in the 
stall ing phase, as was identifi ed in iBaumgardt et al.l (2006) 
and Matsubava shi et ail (12007b . However, the origin of this 
increase is not SDAF (which does not occur for non-migrating 
orbits), but is instead due to the conventional dynamical 
friction being reduced at the periastron (inside the cavity) 
as compared to that at the apoastron (outside the cavity or 
closer to its edge) . Thi s effect was clearly explained in 
iMatsub avashi et al. (2007) and for IMBHs it plays decisive 
role. SDAF does however play an important role if the cusp 
is co- or counter-rotating relative to the IMBH's orbit. 

Observational predictions of SDAF include higher mass ob- 
jects with increasingly larger eccentricities (in non-rotating 
and counter-rotating stellar cusps), and density anisotropy in 
the stellar distribution near the plane of the massive body. 
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